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We report the theoretical discovery of a class of 2D tight-binding models containing nearly-flat 
bands with nonzero Chern numbers. In contrast with previous studies, where nonlocal hoppings 
are usually required, the Hamiltonians of our models only require short-range hopping and have the 
potential to be realized in cold atomic gases. Due to the similarity with 2D continuum Landau levels, 
these topologically nontrivial nearly-flat bands may lead to the realization of fractional anomalous 
quantum Hall states and fractional topological insulators in real materials. Among the models we 
discover, the most interesting and practical one is a square-lattice three-band model which has only 
nearest-neighbor hopping. To understand better the physics underlying the topological flat band 
aspects, we also present the studies of a minimal two-band model on the checkerboard lattice. 
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In fermionic systems, a flat band (a macroscopically 
degenerate manifold of single-particle states) plays an 
important role in the study of strongly correlated phe- 
nomena, due to the vanishing bandwidth. One way to 
achieve a flat band is via the destructive interference 
of electron hoppings, which gives rise to moderately lo- 
calized single-particle eigenstates [IH^I- Landau levels, 
which are formed when a magnetic field is applied to 
a 2D electron gas (2DEG), can be considered as another 
type of flat bands arising in continuum rather than lattice 
2D systems. Different from the examples above, a Lan- 
dau level has a nontrivial topological index (the Chern 
number) . When an integer (or certain fractional) number 
of Landau levels is filled, the system turns into an insu- 
lator with nontrivial topology, known as the integer (or 
fractional) quantum Hall effect (IQHE or FQHE) 0. 

In addition to 2DEG, efforts have also been made to 
realize IQH and FQH effects in lattice systems without 
magnetic field. The first and most celebrated example 
is the anomalous quantum Hall state proposed by Hal- 
dane . More recently, a new class of topological states, 
the time-reversal invariant topological insulators charac- 
terized by a Z2 topological index is discovered in various 
lattice systems (see the recent reviews of Refs. [l^, Il3| 



similar to the FQHE. Recently, there have been some at- 
tempts to find completely fiat bands with nonzero Chern 
numbers in 2D lattice models 14,3. However, it turned 



and references therein). These lattice topological states 
share strong similarities with the IQH states. However, 
the lattice counterpart of the fractional quantum Hall 
states has not yet been discovered. One of the key chal- 
lenges to reach a fractional topological state in lattice 
models lies in the fact that the bandwidth of a topologi- 
cally nontrivial band in these models is usually compara- 
ble or even larger than the band gap. Thus, at fractional 
filling, the system is expected to be a Fermi liquid while 
interaction effects are just subleading corrections. There- 
fore, a flat band with nontrivial topology is expected to 
be the key in realizing lattice fractional topological states 



out that, in all examples except in the quasi-one dimen- 
sional thin torus 14], fiat bands have zero Chern number. 

Since a topological index remains invariant under adi- 
abatic deformations as long as the gap is preserved, a 
straightforward way to form such a flat band is to use 
the spectral flattening trick, i.e., an adiabatic transfor- 
mation from the original Hamiltonian to a new one with 
completely flat bands. This technique is used in the 
classification of topological insulators and superconduc- 
tors 

[BIS- 

However, such a procedure may results 
in long-range hopping making the Hamiltonian nonlocal 
(see Supplementary Material for details). 

In real materials, the exact flatness for a band is 
not a physical requirement and we can relax the con- 
straint a little bit by allowing the band to have a nonzero 
bandwidth but requiring the bandwidth to remain much 
smaller than the band gap. Unfortunately, to the best 
of our knowledge, even such models have never been re- 
ported. In this letter, we propose a generic scheme to 
produce such kind of models based on a special class 
of tight-binding Hamiltonians with short-range hoppings. 
The band structure of these models contains nontrivial 
band touchings with quadratic dispersions (in contrast to 
the linear ones near a Dirac point), which are protected 
by the time-reversal and lattice point-gr oup symmetries 
as well as the nontrivial topology [19|, |20]- When the 
time-reversal symmetry is broken, a band gap opens up at 
the band touching point and the bands can have nonzero 
Chern numbers. By slightly tuning the short-range hop- 
ping strength, we find that some of the topologically non- 
trivial bands can become nearly-fiat. We believe that 
this mechanism is very general and applies to any tight- 
binding model with quadratic band touchings. Surpris- 
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ingly, in some of these models, this nearly-flat band sit- 
uation is found even with only nearest-neighbor (NN) 
hopping. These nearly-flat bands have a strong analogy 
to the Landau levels in 2DEG and thus may set the stage 
for exploring new fractional topological states. 

Topological flat band with extremely short-range hop- 



ping: the square-lattice model — Consider a square lat- 
tice with two space-inversion odd and one space- inversion 
even orbitals per site, e.g. the p^, Py and dj.2_y2 orbitals. 
This model has been demonstrated in optical lattice sys- 
tems [2^. In fc-space, the Hamiltonian is 
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where dj:, p^ j: and Py ^ are the fermion annihilation op- 
erators at momentum k and the lattice constant is set 
to unity. The NN hoppings between various orbitals are 
described by the hopping amplitudes tdd, tpd, tpp and t'^^. 
The 5 term measures the energy difference between p and 
d orbitals. The term (A) breaks the symmetry between 
the states with angular momentum ±1 [px ± ipy). This 
term breaks the time-reversal symmetry and allows the 
Chern number to take a nontrivial value. 

At A = 0, the time-reversal symmetry is preserved 
and two of the three bands cross at the center (corner) 
of the Brillouin zone (BZ). For A > 0, the bands be- 
come gapped. In order to ensure the "flatness" of the top 
band, we require the energies are equal at the V point, 
M point and X point, which implies 5 = —4:tdd + 2tpp -t- 
A - 2tppA/{Atpp + A) and t'^^ = tppA/{4tpp + A). For 
simplicity, we set tdd — tpd = ^pp = 1- By varying A, we 
found that the ratio of bandwidth/band gap is minimized 
(~ 1/20) at A = 2.8. Here the top and the bottom bands 
carry opposite Chern numbers ±1 while the middle band 
has a trivial Chern number. The band structure of this 







model is shown in Figs. 1(b) and 1(c) where the former is 



computed for periodic boundary conditions (on a torus) 
and the latter on a cylinder with two open edges. The 



edge states appearing in Fig. 1(c) confirm the nontrivial 
Chern numbers of the system. 

A two-band model on a checkerboard lattice — The 
model discussed above has three bands. Here we present 
another model with only two bands. A two-band model 
has the following advantages: (1) its band structure is 
much easier to compute analytically; and (2) the Hilbert 
space is much smaller than models with more bands and 
thus numerical studies become easier. However, here we 
need to allow the next-nearest-neighbor (NNN) and next- 
next-nearest-neighbor (NNNN) hoppings. We emphasize 
that single-band models can only have trivial Chern num- 
bers and thus, a two-band model is the minimal model 
to have topologically nontrivial bands. 

Consider a checkerboard lattice with N N (t), NNN {t[, 
t'2), and NNNN (t") hoppings [Fig. |2(a)| . Here, we allow 
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FIG. 1: (Color online) Chiral quasi-flat band in the three- 
band model on a square lattice. Figure (a) shows the lattice 
structure, where each lattice site contains three orbitals and 
the arrows represent the breaking of the time-reversal sym- 
metry. By putting the system on a torus and a cylinder, the 
single-particle energy spectra are shown in Figs, (b) and (c). 
In Fig. (c), chiral edges states (thick lines) are observed. 

the NN hopping to carry nonzero complex phase (i?!'), 
whose signs are shown by the arrows in Fig. 2(a) These 
complex hoppings break the time-reversal symmetry at 
(j) 7^ nn (n e Z) . The Hamiltonian of this model is 

H = -tJ2 e'^-^ (clc, + h.c.) - <Mc^ + ^-c-) 
{«j-> {{id)) 

-t" E (4c, + /i.e.) (1) 

{({«J»> 

where Ci (c|) is the fermion annihilation (creation) oper- 
ator at site i. The NN, NNN and NNNN sites are rep- 
resented by (?,j), {{i,j)) and (((i,j))). The phase factor 
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FIG. 2: (Color online) Chiral quasi-flat band on the checker- 
board lattice. The lattice structure is shown in Fig. (a), with 
arrows and (solid and dashed) lines representing the NN and 
NNN hoppings respectively. The direction of the arrow shows 
the sign of the phase in the NN hopping terms. Two of the 
NNNN hoppings are shown as the dashed curve. Other con- 
vention are the same as in Fig. [T] 

in the NN hopping terms is (jiij = ztcj) with the sign de- 
termined by the direction of the arrows. The hopping 
strength between NNN sites i-^ takes the value of t[ [t'2) 
if the two sites are connected by a solid (dashed) line. 

The checkerboard lattice has two sublattices, and thus 
the Hamiltonian can be written in the momentum space 
as H = — X]fe V'^^V'fci where = {a^, bj^) is a two com- 
ponent spinor and H is a 2 x 2 matrix 



H = [{t[ + i2)(cos kx + cos ky) + At" cos kx cos ky\I 
k k 

-I- 4t cos (bicOS — cos —)(Tx ^ 

+ ^ t'2) {cos kx — cos ky)a. 



k k k k 

4tcos0(cos cos -^)'^x + 4:tsmip{sm sin -^)<^v 

(2) 



Here / and (yx,y,z are the identity and Pauli matrices. 

At = or nTT, the time-reversal symmetry is pre- 
served and the two energy bands cross at the corner of 
the BZ 0. At ^ nTT, a gap opens up between these 
two bands. Due to the breaking of the time-reversal sym- 
metry, each band can carry a nonzero Chern number. In 
order to reach a flat band, we require the energies of the 
top band are equal at at the F point, M point, K point 
and at fc = (±7r/2, ±7r/2). With t = 1 and = 7r/4, 
this condition implies t — \, t\ = —t^ = 1/(2 + \/2), 
t" = 1/(2 2a/2). With these values, the top band be- 
comes very flat, with bandwidth of about 1/30 of the gap 
[Fig. 2(b)| and each of the two bands now carries Chern 
number ±1. We further verify this conclusion via the 
study of the chiral edge mode [Fig. 2(c)| . 



Discussion — In addition to the models discussed 
above, similar effects can be observed in other models 
with quadratic touching. For example, if we allow NN 
and NNN hoppings, both the kagome lattice and the hon- 
eycomb lattice with the Px and Py orbtials [1, [21], [l^] can 
support this type of nearly-flat bands when the time- 
reversal symmetry is broken. Here the kagome-lattice 
model is a three-band one, while the other has four bands. 

Here, we compare the models with quadratic band 
touching and those with Dirac points. Due to fermion 
doubling, the Dirac points need to appear in pairs. In 
order to reach an insulating phase starting from a semi- 
metal with two Dirac points, a nonzero mass need to be 
introduced at each of the two Dirac points. However, 
depending on the relative sign of the two Dirac masses, 
the resulting insulator can be either topologically trivial 
or nontrivial 
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Due to this uncertainty on the topo- 
logical structure, the nearly-flat band from a model with 
Dirac points may be topologically trivial and thus irrel- 
evant to our interests. On the contrary, for the models 
with quadratic band touching, the constraint of fermion 
doubling is absent and there is a single crossing point. 
This crossing point can be regarded as two Dirac points 
merging together, and is protected by time reversal sym- 
metry and discrete rotational symmetry, e.g., C4 in the 
checkerboard lattice model. Since two hidden Dirac cones 
have the same chirality [l^, the energy bands will have 
nontrivial topological numbers once the gap is opened by 
breaking time reversal symmetry by complex hoppings 
which do not break discrete rotational symmetry [23|. 
Therefore, we can focus on the flatness of the band, with- 
out worrying about finding a topologically-trivial band. 
This discussion is also the reason that we considered the 
(i-orbital in the three-band square lattice model. Without 
the d-orbital, the two p-band have two quadratic band 
crossings at the center and the corner of the BZ, instead 
of just one. Therefore, the resulting insulator may be 
topologically trivial. 

Although the band gap is much lager than bandwidth 
in these models, it is not clear whether such nearly-flat 
bands are equivalent to Landau levels in 2DEG. How- 
ever, the Berry curvature (Fig. [3]) in momentum space 
shows no sharp features and the only length scale is the 
lattice constant, in sharp contrast to the cases in which 
the Berry curvature has delta-function like peaks, e.g. 
Ref. [i^j ■ Thus, we argue that the topological nearly- flat 
bands we propose are very similar to 2D Landau levels 
and we expect FQHE at fractional fillings when repul- 
sive interactions are turned on. We note in this context 
that even 2D Landau levels have a short lattice length, 
typically 10-100 times smaller than the magnetic length, 
underlying the real physical 2DEG. A more detailed nu- 
merical study will be presented in our future work. 

When spin degrees of freedom are taken into account, 
the discussion above can be generalized to the time- 
reversal invariant Z2 topological index by just substi- 
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FIG. 3: (Color online) Distributions of the Berry curvature 
in momentum space for the flat bands in (a) the square-lattice 
model and (b) the checkerboard lattice model. 

tuting the time-reversal symmetry breaking terms into 
corresponding spin-orbit couplings (where spin up and 
down particles break the time-reversal symmetry in the 
opposite way and thus the time-reversal symmetry is re- 
covered when both spin species are taken into account). 
In such a way, it may even be possible to realize fractional 
topological insulators 25| in these models. 

Experiment realization — In the discussion above, we 
only provided the optimum values for the parameters at 
which the flatness of the band is maximized. However, 
the nearly-flat band does not require strictly fine-tuning 
to reach. In fact, even if the parameters are changed 
by about 10%, the band remains to be fairly flat in the 
models we studied (see Supplementary Material for de- 
tails). Because of this stability and the simpleness of 
these models, we believe that experimental realizations 
of these models are possible in both condensed matter 
systems and ultra-cold atomic gases. The insulating gap 
in these systems can be opened via spontaneous symme- 
try breaking if a small amount of short-range repulsions 
are introduced 1^ 2^, 2^. The same effect can be ex- 
pected via explicit symmetry breaking, e.g. by introduc- 
ing a magnetic field (for charged particles) or an artificial 
gauge field (for charge neutral particles) well as 

by rotating the lattice 22]. In recent experiments, some 
optical lattices have been constructed whose band struc- 
tures are described by the square-lattice model a nd the 
honeycomb-lattice model discussed above 
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sidering the fact that hopping strength can be tuned rela- 
tively easily in cold gases via varying the optical lattices, 
these cold-atom systems may be the leading candidates 
for the realization of the topological physics predicted in 
our work. 
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Note added. — Related work has recently been done in 



Refs. [32, |33| . Very recently, the existence of fractional 
quantum Hall effect in our model has been confirmed by 
exact numerical studies as reported in Ref. 34 1. 
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I. TOPOLOGICAL FLAT BAND FROM LONG-RANGE HOPPING 

If nonlocal terms are allowed, flat band with nonzero Chern number can be constructed using an idea similar to 
the spectral flattening trick. For an n-band model whose bands carry nontrivial Chern numbers, e.g., Haldane's 
honeycomb lattice model, we can write the Hamiltonian in momentum space as an n x n matrix 

n{k) = U{k)A{k)U\k), (SI) 

where A(fc) is the diagonal matrix with Ea{k) (a = I, n) on the diagonal and U{k) is a unitary matrix. From A(fc), 
we can construct a new diagonal matrix A(fc) in which the energy of the topologically non-trivial band, say Eb{k), is 
replaced with Eb, the constant independent of k. Then we construct a new Hamiltonian: T-L{k) = U{k)A{k)W (k). It 
is obvious that 'H(fc) has a completely flat band at E = E^, which carries nonzero Chern number since the eigenvector 
corresponding to Eb remains unchanged from that of Eb{k). Note that as long as the 6th band is isolated from the 
other bands, the other energy levels in A(fc) can also be modified. One can obtain the Hamiltonian in real space 
by the inverse Fourier transform of 'H{k). However, in this construction, the hopping amplitudes in general remain 
nonzero even between sites with arbitrary large separation. 

II. FLATNESS OF THE BAND 




/pd t'/t 

(a) (b) 

FIG. SI: (Color online) The ratio of band gap/bandwidth for the nearly-flat band in (a) the square-lattice model and (b) the 
checker-board lattice model. In Fig. (a), the horizontal and vertical axes are tpd and A, correspondingly. For other parameters, 
we set 5 = —Udd + 2tpp + A — 2tppA/{4:tpp + A), t'pp = tpp/{Atpp + A) and tdd = tpp = 1- In Fig. (b), we set t[ = = t' and 
(j> = 7r/4 with the two axes being t'/t and t" /t. The dot at the center of each figure marks the parameters used in the Letter. 
In both figures, the parameters are varied by the amount of ±20%. The ratio of band gap/bandwidth are marked in the figures 
for each contour. 

In Fig. ISH we show the contour plot of the ratio of band gap/bandwidth for the nearly-flat band at different 
parameters (±20% away from the values we used in the Letter). As can be seen from the figures, in a wide region of 
the parameter space, this ratio remains large indicating a very flat band. 

We'd also like to emphasize that in this Letter, we do not intend to maximize the flatness of the band, but to 
demonstrate a generic technique which provides nearly-flat topological bands. For example, in the square-lattice 
model, we set tdd = tpd = tpp = 1 in the Letter for simplicity and maximized the ratio of band gap/bandwidth 
within this constrain by varying A, which results in a very-flat band with nontrivial topology. However, as shown in 
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Fig. ISll fa). it is possible to enhance this ratio by removing the requirement of tdd = tpd = tpp = 1, but it is not the 
central focus of this work to find the optimum value for each control parameter to maximize this ratio. 



